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THE  INFLUENCE  OF  QUANTIZING  AND  SAMPLING 
ON  THE  SIGNAL-TO-NOISE  RATIO 
OF  A  COMPLEX  CORRELATOR 


By 

T.  Kooij  and  R.  Laval 

ABSTRACT 


The  combined  effects  of  sampling  and  phase-quantizing  of  phase- 
modulated  signals  in  Gaussian  noise  are  studied  and  expressed  in 
terms  of  the  output  signal-to-noise  ratio  of  a  matched  filter 
correlator.  The  obtained  results  are  the  basis  for  finding  the 
sampling  and  quantizing  losses  in  other  types  of  correlators 
or  filters.  Since  only  the  signal  bandwidth,  and  not  the 
carrier  frequency,  is  considered  to  be  of  importance,  complex 
envelopes  are  used  to  describe  the  signals,  leading  to  complex 
correlator  types.  The  results  are  summarized  in  three-dimensional 
graphs . 
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LIST  OF  SYMBOLS 


a 


attenuation 


A(t) 

B(t) 

c(t) 


E 


cosine  component  of  expected  signal 

sine  component  of  expected  signal 

correlation  function  of  cosine  components  of 
received  and  expected  signals 

expectation,  statistical  average 


F 


central  frequency  of  narrow-band  signal  before 
demodulating 


weighting  function  of  complex  matched  filter 

average  length  of  normalized  noisy  input  vector 
for  given  aS  and  N(t) . 

probability  that  clipped  component  of  input  signal 
has  the  same  sign  as  the  same  component  of  the 
expected  signal 


,  •  sin  tt 

i  2 - E_ 

I*  -CD  P 


.  —  .  arcsin 
TT 


P 

sin  ^  tt 


noise  factor 


m 


m 


m  sampling  rate  multiplier  (m  =  1  for  sampling  with 

frequency  AF),  or 

phase  quantizing  factor  (m  =  2  for'^-  phase  quantizing) 

* 


2 


n(t) 

N 

N(t) 

N 

p[«Pi»  <Pa] 
Rs(t) 

S 

SNR 

T 

a(t) 

P(t) 
y(  t) 

r(  u>) 

6(t) 

AF 

At 


input  noise  signal  per  component 
T  .  AF  compression  factor;  gain  factor 
modulus  of  complex  input  noise  signal 
average  noise  modulus 

joint  probability  density  function  of  cpn ( t x  ) 
and  cpn  ( t a  ) 

time-autocorrelation  function  of  one  component 
of  the  expected  signal 

constant  modulus  of  expected  signal 

signal- to-noise  ratio 

time  duration  of  expected  signal 

cosine  component  of  received  signal 
sine  component  of  received  signal 
complex  correlation  function 

Fourier  spectrum  of  complex  correlator  output 

Dirac  function 

signal  and  noise  bandwidths 

sampling  interval  =  — 

H  m .  AF 
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tl  a 


T|(t) 

H(w) 

V(t) 


U  t) 


expected  signal  =  A(t)  +  j  B(t)  =  S  e 
Fourier  transform  of  expected  signal 

jcpn(t) 

complex  input  noise  signal  =  N(t)  .  e 

j[«ps(t)+cp0] 

received  signal  =  a( t )+j p ( t )  = aS . e  +  v(t) 

Fourier  transform  of  input  noise 
phase-quantizing  step  (m  £  2) 


T 


CPo 

<us(t) 


Vt) 


autocorrelation  function  of  the  input  noise 


power  per  expected  signal  component  =  £sa 


power  per  input  noise  component 


correlation  delay  parameter 


phase  difference  between  received' and  expected  signal 


phase  angle  of  expected  signal 


phase  angle  of  complex-input  noise 
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INTRODUCTION 


Correlators  belong  to  the  class  of  devices  that  average  the 
product  of  two  functions.  Although  for  a  correlator  these  functions 
are  signals  in  time,  they  have  the  operations  of  product-averaging 
in  common  with  other  "signal  processors"  that  mathematically  perform 
the  same  averaging  over  a  different  variable,  for  example  over 
azimuth  in  the  case  of  a  bearing  estimator.  All  systems  deriving 
decisions  from  second-order  statistics  perform  basically  the  same 
operations  as  a  correlator. 

The  name  "correlator"  thus  defined  varies  with  the  properties  of  the 
functions  over  which  the  product  is  averaged.  If  both  functions 
are  completely  deterministic,  the  device  is  only  a  calculator, 
calculating  the  mean  product.  For  non-triviality  we  shall  require 
therefore  that  at  least  one  of  the  two  signals  is  disturbed  by  some 
non-deterministic  corruption.  If  one  signal  is  corrupted  and  the 
other  not,  the  latter  is  called  the  "filter  function"  or  "system 
function",  or  in  special  circumstances  the  "weighting  function"  or 
"impulse  response  function" .  The  properties  of  this  second  signal 
define  the  name  of  the  system,  i.e.  "matched  filter"  if  the  second 
function  resembles  the  non-corrupted  input  signal,  simply  "filter" 
in  the  case  of  an  unspecified  function,  or  "frequency  analyser" 
if  the  second  function  is  a  pure  sine  wave. 
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Only  if  both  of  the  averaged  signals  are  disturbed  by  a  random 
background  does  the  averaging  device  properly  deserve  the  name 
"correlator" .  The  table  below  may  help  to  distinguish  between 
the  possible  correlator  systems. 


Type 

of 

Correlator 

Signal  x(t) 

Signal  y(t) 

Averaged  of  shifted 
time  product 
x(t) y(t-T)  dt 

l 

Known 

Random 

Known 

Random 

i 

X 

X 

calculator 

2 

X 

X 

Filter 

3 

X 

known 
part  of 
x(t) 

matched  filter 

4 

X 

j  2nft 

e 

f requency  analyzer 

5 

X 

X 

proper  correlator 

All  correlators  have  in  common  the  need  for  a  memory  for  both 
signals  involved.  In  an  analogue  filter  one  of  the  signals  is 
stored  in  the  configuration  of  the  filter  components,  the  other 
signal  (the  input)  has  its  running  past  stored  as  partial 
contributions  to  the  electromagnetic  fields  in  the  components. 

Modern  signal-processing  has  for  the  last  ten  years  made  extensive 
use  of  binary  digit  memories  (cores,  delay  lines,  flip-flops). 

The  cost  of  these  memories  is  generally  presented  per  bit. 
Reduction  of  the  necessary  number  of  bits  therefore  is  and  will  be 
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of  increasing  economic  importance.  Digitizing  means  sampling  in 
time  and  quantizing  in  amplitude,  and  this  paper  tries  to  present 
a  systematic  approach  to  describe  the  effects  of  economizing  on 
both  of  these  operations  when  applied  to  phase-modulated  signals 
corrupted  by  Gaussian  noise,  thereby  permitting  the  cost  per  bit 
to  be  weighed  against  the  loss  of  performance. 

Rather  than  studying  the  effects  of  bit-reduction  on  all  systems 
given  in  the  table  above,  it  has  been  decided  to  study  the  effects 
on  only  one  signal  at  a  time.  This  is  done  by  comparing  or 
correlating  the  quantized  signal  with  its  uncorrupted  replica. 

The  results  of  this  study  are  therefore  only  directly  valid  for 
a  matched  filter,  and  only  valid  by  implication  for  the  other 
correlator-types . 

Another  restriction  in  this  study  is  the  requirement  of  pure 
phase-modulation  of  the  signals.  This  makes  the  results  less  general 
in  theory,  but  in  practice  many  signals  applied  in  systems  for 
active  detection  and  acquisition  (e.g.  tracking)  are  of  this 
constant-amplitude,  phase-modulated  type. 

The  expression  "complex"  .  in  the  title  indicates  that  each  signal 
has  been  considered  to  be  demodulated  by  its  central  frequency, 
leaving  a  "complex  envelope"  (in  practice  a  pair  of  in-phase  and 
quadrature  components) .  The  correlator  for  such  signals  naturally 
bears  the  name  "complex  correlator" .  If  the  original  phase- 
modulated  signal  is  demodulated  with  a  frequency  different  from 
the  carrier,  the  complex  envelope  will  contain  a  residual  modulation 
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equal  to  the  difference  frequency.  Because  of  the  uniform  phase 
distribution  of  the  modulating  inter-frequency,  the  signal-to-noise 
ratios  derived  in  this  pape.r  for  the  purely  complex  correlator 
apply  equally  well  to  a  correlation  system  that  for  technical 
reasons  uses  a  residual  carrier  (ac  correlator) . 

A  widely  used  method  of  amplitude-quantizing  is  infinite  clipping, 
which  destroys  all  amplitude  information  and  retains  only  the  real 
axis  crossings  of  the  phase.  By  linear  interpolation  between  the 
crossings,  a  phase  curve  may  be  reconstructed  that  will  resemble 
the  phase  of  the  original,  non-clipped  signal.  Although  there  are 
different  approaches  of  clipped  signal  reconstruction  (e.g.  using 
Fourier  transform  methods) ,  in  this  paper  the  theoretical 
reconstruction  method  is  to  interpolate  the  phase  between  the  real 
axis  crossings,  as  indicated  above,  and  to  consider  this  function  to 
be  the  phase  of  a  signal  with  a  constant  modulus,  so  that  the 
reconstructed  signal  is  purely  phase-modulated. 

It  will  be  clear  that  this  method  is  much  more  unfavourable  to 
amplitude-modulated  signals  (which  are  completely  stripped  of  their 
side-bands  by  the  clipping)  than  to  phase-modulated  signals. 

Gaussian  noise  has  its  information  more  or  less  equally  divided 
between  pure  phase  and  amplitude,  so  that  phase-modulated  signals 
plus  noise  are  in  the  favoured  half. 

If  the  clipping  and  interpolation  process  is  thought  to  have 

happened  at  a  sufficiently  narrow  relative  bandwidth,  the  reconstructed 

wave  will  have  a  phase  curve  that  closely  follows  the  original  phase 
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law,  so  that  the  same  will  be  true  for  the  phase  of  the  complex 
envelope  after  demodulation  with  the  carrier.  The  total  effect  will 
be  that  the  complex  low  frequency  signal  has  been  trimmed  or 
"shaved"  to  a  constant,  normalized  modulus. 

This  is  the  only  method  of  modulus  quantizing  that  is  considered 
in  this  report.  It  is  clear  that  such  a  system  is  not  able  to 
provide  information  about  the  absolute  energy  of  the  input  signal, 
but  only  about  energy  ratios  with  respect  to  other  inputs,  such  as 
noise . 

The  other  type  of  quantizing  considered  is  phase-quantizing  of 
the  complex  envelope.  Although  phase-quantizing  is  equivalent  to 
amplitude-quantizing  of  the  in-phase  and  quadrature  components, 
these  quantizations  are  not  independent,  but  related  to  each  other 
as  cosine  and  sine. 

It  is  known  that  a  proper  Nyquist  sampling  frequency  (twice  the 
signal  bandwidth)  is  sufficient  to  describe  a  signal  without  loss 
of  information.  However,  if  this  same  sampling  frequency  is 
tentatively  used  on  the  same  signal  after  clipping,  the  losses  — 
due  to  mismatched  sampling  frequency  as  a  result  of  the  widened 

t 

signal  bandwidth  —  will  generally  be  even  greater  than  they  already 
were  due  to  modulus  clipping  alone. 

After  an  explicit  summary  of  constraints  and  assumptions,  the  effects 
on  correlator  performance  of  combined  modulus  normalizing  (clipping 
the  high  frequency  input) ,  phase-quantizing,  and  sampling,  will  be 
followed  step  by  step. 

A  summary  of  the  results  is  given  in  graphical  form  in  the  last 
chapter . 
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1. 


ASSUMPTIONS  AND  DEFINITIONS 


The  complex  low-frequency  signals  considered  in  this  report  will  be 
assumed  to  have  the  following  properties: 

a.  The  time  duration  of  the  signal  is  T  (seconds). 

b.  The  signal  has  a  flat  spectrum  with  bandwidth  AF(Hz), 
centred  at  zero  frequency. 

c.  The  signal  has  a  uniform  phase  distribution. 

d.  The  expected  signal  has  a  constant  modulus  S  =  G  */"~2  . 

e.  The  received  signal  is  an  attenuated  copy  of  the  expected 
signal,  with  a  constant  modulus  as  and  with  a  constant 
phase  shift  qp0  with  respect  to  the  expected  signal . 

f.  The  expected  signal  is  not  quantized. 

g.  The  quadrature  components  of  the  noise  have  a  flat 

spectrum  (bandwidth  AF,  centred  at  zero  frequency)  and 

a  Gaussian  amplitude  distribution  (variance  a3  ), 

n 

h.  The  input  SNR  will  be  defined  to  be  the  average  signal 
power  divided  by  the  average  input  noise  power. 

i.  The  output  SNR  will  be  defined  to  be  the  square  of  the 
average  magnitude  of  the  correlation  peak  divided  by  the 
average  output  noise  power. 
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2. 


QUANTIZING  AND  SAMPLING 


2 . 1  Quantizing 

The  different  types  of  quantizing  of  the  received  signal  considered 
♦ 

in  this  report  are: 

a.  Normalizing  the  modulus  to  a  constant  length  with 
conservation  of  phase  information. 

b.  Phase-quantizing  combined  with  normalizing  the  modulus. 

This  operation  converts  the  signal  to  a  vector  in  the  complex 
plane  with  a  fixed  magnitude  and  a  finite  number  of  possible  phase 
angles.  These  angles  will  usually  be  chosen  regularly,  so  that 

a  regular  vector  star  may  be  Used  to  describe  this  kind  of 
quantizing. 

The  simplest  form  of  phase  quantizing  is  a  vector  star  consisting 

If  the  star  is  orientated  at  45° 
to  the  horizontal,  this  form  is 
equivalent  to  clipping  the 
quadrature  components,  as  may  be 
done  in  a  technical  system  (Ref.  1  )  . 
The  limiting  case  of  infinite  phase- 
quantizing  makes  type  (b)  euqal  to 
type  (a) :  complete  phase 
description  and  normalized  modulus. 


of  four  vectors  (Fig,  1). 
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2 . 2  Sampling 


The  sampling  theorem  states  (Ref.  2)  that  a  high-frequency  signal 
with  bandwidth  AF  is  described  completely  by  double-samples 
taken  at  time  intervals  of  l/AF  seconds.  Each  sample  determines 
the  amplitude  of  the  signal  and  its  phase  with  respect  to  the 
carrier,  and  could  therefore  be  represented  as  a  complex  sample. 

The  sampling  frequency  AF  is  also  found  by  considering  that 
demodulating  the  high-frequency  signal  translates  the  spectrum  to 
zero  frequency  (Fig.  2),  The  two,  real,  low-frequency  signals 
describing  the  envelope  and  the  phase  of  the  modulating  function 
will  have  a  physical  bandwidth  of  ^AF,  thus  requiring  a  sampling 
frequency  of  AF  Hz,  as  before. 


AF 


I - - 1  Demod. 

I  !  I  - - 

J _ I  I _ 

-F 


Spectral 

density 


AF  O  _AF 
2.  2 

Fig.  2 


Demod. 


AF 


4- 

F 


frequen  cy 


Quantizing  operations  —  such  as  clipping  the  low-frequency 
quadrature  components,  which  in  itself  causes  loss  of  information 
and  deterioration  of  the  SNR  in  detection  —  tend  to  broaden  the 
bandwidth  of  the  original  signals.  The  sampling  frequency  AF 
therefore  will  no  longer  be  adequate  to  describe  the  quantized 
signals  completely,  resulting  in  even  more  losses  in  the  SNR  of  a 
detector . 
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It  is  to  be  expected  that  the  losses  will  diminish  if  the  sampling 
frequency  is  increased.  This  is  seen  clearly  in  the  case  of  the 
clipping  of  the  quadrature  components,  which  are  then,  after 
clipping,  described  completely  by  their  zero-crossings.  These 
zero-crossings  can  be  determined  with  a  higher  accuracy  when  the 
sampling  frequency  is  increased.  The  losses  due  to  the  clipping 
alone  are  found  by  increasing  the  sampling  rate  to  continuous 
description  in  time  ('infinite  sampling'). 
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3  .  THE  REFERENCE  SYSTEM 

The  behaviour  of  a  complex  correlator  without  any  quantizing  or 
sampling  will  be  taken  as  a  standard  of  performance. 

With  the  assumptionsof  Chapter  1,  and  using  Greek  symbols  for  complex 
functions,  we  consider  the  complex  correlation 


y(t)  =  l/T 


£(t)  T)  (t  -  t)  dt 


(Eq.  1) 


where 


v(t)  =  received  signal 


expected  signal 


V(t)  -  N(t)  eiCpn(t) 


noise 


3  .  1  Signal  Correlation  Peak 


This  correlation  function  yields  for  T  =  0  (coincidence): 


y(0)  =  y  (0)  +  Y  (0) 

s  in 


(Eq.  2) 


(aS2/T)eitp»  rei[‘»s<t)-'Mt)llt+s/T  dt 


u 
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The  actual  value  of  y(0)  depends  on  the  contribution  y^(0)  due  k° 
the  noise.  The  mean  value  of  y(0)  is  found  by  taking  the 
statistical  average 


E  [v(0) ]  -  ys(0)  +  E  [Yn(0)J  , 


(Eq.  3) 


where  y^O)  is  the  cmly  random  variable,  and 


(Eq.  4) 


The  integrand  is  zero,  being  the  average  of  a  vector  with  uniform 
phase  distribution.  Therefore 


E  [Yn(0>]  =  0 


and 


(Eq.  5) 


It  is  seen  that  the  initial  phase  difference  cp0  between  received 
and  expected  signal  is  repeated  in  the  phase  of  the  correlation  peak. 

3 . 2  Noise 

The  average  power  of  the  correlator  output  if  the  input  is  noise  only, 
is  found  by  considering  the  complex  correlator  as  a  matched  filter 
i.e.  a  filter  having  a  weighting  function  h(t)  equal  to  the  complex 
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conjugate  of  the  expected  signal,  running  backwards  in  time: 


h(t)  =  (i/T)  n*(-t) 


(Eq.  6) 


The  factor  ( l/T)  is  due  to  the  normalizing  factor  in  the  correlation 
( Eq •  1 )  • 


The  filter  frequency  response  is  equal  to  the  complex  conjugate  of 
the  spectrum  of  the  expected  signal .  Written  as  a  Fourier  pair 


(l/T)  n*(-t)  ^  (  l/T)  H*(U)) 


(Eq.  7) 


if  T)  (t)I^H(tD) 


(Ref. 2  ,  rules  2  &  3) 


It  is  further  known  (Ref.  3  )  that  the  spectral  density  of  the 
output  of  a  linear  filter  is  equal  to  the  input  spectral  density 
multiplied  by  the  squared  modulus  of  the  filter  response: 


|rn(u>)|  -  | 


Z  <*)  I  (  l/T3)  |H*(u)  I 


where 


|T  (U>)  |  =  output  noise  spectral  density 


input  noise  spectral  density 


i  h*(®> 

T 


«*  matched  filter  response 


(Eq.  8) 
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The  total  input  power  is  equal  to  20 *  (two  orthogonal  noise 

signals  of  power  o ,  causing  the  modulus  to  have  a  Rayleigh 

distribution  with  variance  2o 3  (Ref.  3  ).  This  input  power  is 

n 

uniformly  spread  out  over  AF ,  so  that  the  input  density  spectrum 
has  a  height  equal  to  20^  / AF  (Fig.  3)* 


Input  Filter  Ou.tpu.b 

Fig.  3 


■ft  * 

The  filter  factor  (l/T2)  |  H  ( U))  |  in  Eq.  8  is  proportional  to 

2 

the  energy  density  spectrum  |H(u))|  of  the  expected  signal,  having 

a  total  energy  equal  to  S2T  =  2o2T,  also  uniformly  spread  out 

s 

over  the  bandwidth  AF .  The  height  of  the  squared  filter  response 
is  thus  (see  Fig.  3): 


(l/AF)(l/T2)  |  H'^(UJ)  |  =  (  l/Ta)  (2o  2  T)/AF  =  20  2  /TAF  =  20  2  /N 

s  s  s 

(Eq.  9) 

N  =  T  AF  =  compression  or  gain  factor. 
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3  N> 


The  total  output  power  is  found  by  integrating  the  output  density 
spectrum  over  AF,  yielding  (see  Fig.  3)5 


0 


3 

n( out ) 


3 

IT  (  U>)  I  AF  =  (  2 a  3 / AF )  (2o  3  /N)  AF  = 
1  n  1  n  s 


4c  3  a  3 /N 
s  n 


(Eq.  10) 


3 • 3  Signal-to-Noise  Ratio 


The  output  SNR  of  the  correlator,  being  the  ratio  between  the 
squared  average  of  the  correlation  peak  ( Eq .  5)  and  the  mean  output 
noise  power  ( Eq .  10)^  becomes 


(SNR) 


ref 


■Uaac»)  /(4o»  a  2/N)  = 

D  oil 


=N ( SNR)  . 

in 


(Eq.  11) 


It  is  seen  that  the  gain  of  the  system  is  equal  to  the  dispersion 
or  compression  factor  N  =  T . AF . 

If  we  write  the  result  in  Eq .  11  slightly  differently,  we  obtain 

the  well-known  relation 


(SNR)  .  =  ( 2aaa 2 T)/( 2a 2  /AF)  = 
ref  s  n 


total  received  signal  energy 

noise  power  per  Hz 


(Eq.  12) 
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4 


CALCULATION  OF  THE  COMPLEX  SNR  BY  REAL  AND  IMAGINARY  PARTS 


If  the  complex  functions  are  described  by  their  real  and  imaginary 
parts  instead  of  by  the  vector  notation  of  the  foregoing  chapter, 
the  complex  correlation  of  Eq.  1  becomes: 


y(t)  =  l/T 

1 

=  l/T 
where 


£(t)n*(t-T)dt=l/T  J  [<x(t)+jp(t)j  LA(t-T)-jB(t-T)  J  dt  = 


(aA  +  pB)  dt  +  j  l/T  (pA  -  aB)  dt 


f 


(Eq.  13) 


£(t)  =  d(t)  +  j  p(t) 
T|(  t )  =  A(t)  +  j  B(t) 


=  received  signal 


expected  signal 


J 


(Eq.  14) 


If  cp0  ,  (the  phase  difference  between  received  and  expected  signal), 
equals  zero,  a  and  p  will  be  attenuated  copies  of  their  expected 
counterparts  A  and  B,  immersed  in  Gaussian  noise.  Only  the  real 
part  of  the  complex  correlation  contributes  to  the  signal  correlation, 
while  the  imaginary  part  in  Eq .  a  13  has  a  mean  value  of  zero 
because  of  the  orthogonality  of  P,A  and  a,B  if  CPo  =  0  • 


Assuming  that  l/T 


aA  dt  (which  ,  for  cpD  **  0,  is  the  correlation 


of  a  real  function  with  its  replica)  has  a  SNR  equal  to  S/N 
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then  the  SNR  for  l/T 


( CXA  +  gB)  dt  will  be  2s/n*/2  because 


the  signal  strength  is  doubled,  while  the  noise  power  increases 

only  by  3  dB,  since  the  quadrature  components,  a  and  g  .  , 

noise  noise 

are  uncorrelated.  The  same  increase  (3  dB)  in  noise  power  results 
from  the  vertical  or  imaginary  contribution  of  the  complete 
correlation  (Eq.  13),  while  the  signal  part  remains  unaltered, 

Cp0  =  0.  The  total  SNR  thus  becomes  2S^N»/2A/2  =  S/N,  which  i 


if 


is 


equal  to  the  SNR  that  was  assumed  for  l/T  OtA  dt . 


J 


If  cp0  /  0 ,  a  and  g  are  no  longer  replicas  of  A  and  B.  However, 
the  whole  system  is  rotationally  invariant  in  amplitude,  due  to  the 
assumed  uniform  phase  distributions  and  the  complete  phase 
description  of  the  expected  signal.  The  results  found  for  qpc  =  0 
are  therefore  generally  valid  as  far  as  amplitudes  are  concerned. 


EXAMPLE; 

The  output  SNR  of  the  reference  system  will  be  determined  by  the 
method  described  above. 


a .  Signal  Correlation  Peak 
If  cp0  =  0,  the  cosine  component  of  the  input  is 

a(t)  =  a  A(t)  +  n(t)  (Eq.  15) 

a  =  attenuation 

A(t)  =  cosine  component  of  expected  signal 
n(t)  =  cosine  component  of  input  noise 
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Correlation  with  A(t)  gives 


C(t)  =  l/T  |_a  A(t)  +  n(t)  I  A(t  -  t) 


dt 


(Eq.  16) 


In  coincidence  (t  =  0)  the  correlation  peak  is  obtained 


C(0)  =  C  (0)  +  C  (0)  =  a/T 
s  n 


A2  (t )  dt  +  l/T 


n(t)A(t)  dt  (Eq.  17) 


The  mean  value  of  the  peak  is 


E  j_C(0)  J  =  Cg(0)  +  E  [cn(0)] 


(Eq.  18 ) 


Cn(0)  being  the  only  random  variable. 


E  [Cn(0)]  =  l/T  E  [n(t)j  A(t)  dt  =  0 


(Eq.  19) 


The  integrand  in  Eq.  19  equals  zero,  because  the  mean  value 
of  the  noise  is  zero . 

Thus  ,  using  Eq .  17, 


r— 

ac  2  T 

|_C ( 0 )  J  -  Cs(0)  =  a/T 

|  A2 (t )dt  —  ST 

a  a  2 
s 

(Eq.  20) 
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where 


a  2  T  =  signal  energy  per  component 

S 


b .  Noise 

In  the  absence  of  signal,  the  correlator  output  at  a  given  instant 
1 1  equals 


Cn(t,)  =  1/T 


n(t)  A(t)  dt 


(Eq.  21) 


This  integral  is  a  stochastic  variable  with  zero  mean  (Eq.  19). 
Its  variance  is  the  statistical  average  of  its  square: 

E[Cn  ]  =  1/Ar2  J  E[n^^  A(t)  n(0)  A(0)~j  dt  d0  = 


1/TS 


■»  p 

A(t)  A(e)  E|  n(  t )  n(  0)  j  dt  d0 
J  J 


(Eq.  22) 


where 


E^n(t)  n(t)  n(  6)  prob^n^  ,ng  j  dt  d6=pn(t-0)  =  pn(  t)  (Eq.  23) 


in  which  p  (t)  =  autocorrelation  function  of  the  noise. 

Kn 

By  a  change  of  variables  t=t,  0  =  t  -  T,  the  integral  of  Eq.  22 
becomes : 


E(C  *  )  =1/T: 


A(t)A(t-T) Pn ( t) dt  dT  -1/T 


r  P 

dTpn(  T)  j_l/T  A(t)A(t-f)dt  j 


(Eq.  24) 
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The  integral  between  the  brackets  in  Eq.  24  is  the  time-autocorrelation 

function  R  ( t)  of  the  expected  signal  A(t)  .  Due  to  the  flat 

spectrum  of  A(t),  its  autocorrelation  function  is  approximately 

a  sin  x/x  function.  Since  A(t)  =  0  for  |t|>^T,  the  correlation 

function  R  (t)  *=  0  for  |t|>T.  The  limited  time  duration  of  the 
s 

signal  causes  a  tapering  or  R  (  t)  for  larger  values  of  t.  Therefore, 

s 

if  TAF>>1,  we  may  approximate  the  time-autocorrelation  function  by 


Rg(T)  =  1/T  I  A(t)A(t-T)dt  =  0s^l-(  I  T  l/T)  j^(sinnAFT/TTAFT)  j  dT  (Eq.2S) 


Equation  24  now  becomes 

T 


(sinTTAFT/nAFT-)  dT  (Eq.  26) 


Since  the  contribution  to  the  integral  is  small  for  large  values 
of  t  (assuming  that  TAF>>1),  the  tapering  factor  ^l-(|t|/T)j 
will  be  neglected;  so  that 


T 


2 


E  (C  2  )  =  1/T 

n 


-T 


nT  AF 


2 


=  (  1/T  AF)  o  2  o  2 
n  s 


-nT  AF 


(Eq.  27) 
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The  integral  limits  may  now  be  extended  from  minus  to  plus 
infinity,  since  the  tails  of  the  integral  contribute  little  to 
the  result  if  T  AF  >>1. 


The  value  of  the  integral  from  -•  to  +•  is  equal  to  rr,  so  that 


E(C  8)  =  ( 1/T  AF)  o  8  0  3 
n  n  s 


1/N  (o  3  0  3) 
n  s 


(Eq.  28) 


Equation  28  gives  the  statistical  average  of  the  square  of  the 
correlator  output,  if  the  input  is  white  noise  with  power  equal 
to  an3  and  bandwidth  AF.  Because  of  the  property  of  ergodicity, 
this  statistical  average  is  equal  to  the  time  average  or  power  of 
the  output  noise. 

c .  Signal-to-Noise  Ratio 

From  Eqs.  20  and  28  it  xs  possible  to  determine  the  SNR  of  the 
output : 

(SNR)  =  [~aao4  /(  l/N)o  8  0  3  1  =  N(a3o  3  /o  3)  (Eq.  29) 

ALs  nsj  sn 

It  has  been  shown  in  the  beginning  of  this  chapter  that  the  output 
SNR  of  the  complete  complex  correlation  is  equal  to  the  value  found 
in  Eq.  29  for  the  correlation  of  the  real  signal  components  only. 
Although  the  calculation  has  been  made  for  cp0  **  0,  the  result 
is  also  valid  if  cp0^0. 
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Thus 


(SNR) 


ref 


N(aaa  3 /a a) 
s  n 


N. (SNR) . 

m 


See  also  Eq.  11. 


(Eq.  30) 
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5 


THE  SAMP LED  REFERENCE  SYSTEM 


Sampling  the  complex  input  signal  is  equivalent  to  sampling  its 
quadrature  components  synchronously. 

If  the  sampling  interval  is  At,  the  sampled  cosine  component 
(indicated  by  a  dot  on  the  non-sampled  function)  may  be  represented 

by 


&(t)  =  E  d(tk)  At  6(t  -  t^) 


(Eq.  3D 


The  scaling  factor  At  is  introduced  to  give  the  sampled  signal 
the  same  impulse  as  the  non-sampled  signal.  Correlation  with  A(t) 
gives 


c&(t)  =  (  1/T)  J  &(t)A(t  -  T )  dt  =  (  1/T)  E  a(tk)At 


A(t-T)  6(t-t  )  dt 

IV 


=  (1/T)  E  a(t  )  A( t  - T )  At  (Eq.  32) 
k  R  k 


5 . 1  Signal  Correlation  Peak 


If  cp0  =  0  and  if  the  signals  are  in  coincidence  (t  =  0),  the 
correlation  peak  is 


C^(0)  =  (1/T)  E  [a  A(tk)  +  n(tk) j  A(tk)  At 
k 


(Eq.  33) 
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The  mean  value  is 


e[c^(0)  J  =  +  E[c,;(0)}=(a/T)S  A2(tk)  At  +(  l/T)E  E[n(^k)  j  A(tfc)  At 


(Eq.  34) 


The  statistical  average  of  each  noise  sample  E^n(tk) j  equals 

,  so  that  J  =  ^  and 


zero 


Efc.(0)l  =  C*(0)  =  (a/T)  E  A3  ( t  )  At 
L  Gt  J  s  k 


(Eq.  35) 


If  the  number  of  samples  is  large  enough,  the  value  of  E  A2(t  )At 


is  approximately  equal  to 


A3(t)  dt  =  a  3  T  =  energy  of  the 
s 


expected  signal.  The  actual  value  depends  slightly  on  the 
position  of  the  sampling  comb. 


The  correlation  peak  (Eq.  35)  of  the  sampling  system  will  thus 

have  the  value  ao  2  ,  as  did  the  peak  of  the  non-sampling  system, 

s 

(Eq.  20). 


5 . 2  Noise 

The  power  of  the  output  noise  is  found  by  determining  the  variance 
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of  the  correlator  output  if  the  input  is  noise  only. 


E[(>3(tt)j=  (1/T3)  ZZA(tk)  A(ti)  E[n(tk)  n(t^)J  At2  (Eq.  36) 

k  t 

with 

E[n(tk)n(tt)J  =  Pn(tk  -  tt)  =  pn(p  .  At) 

0  (  t)  =  0  s  (sin  TT'AF* t/tt-AF*  t)  =  autocorrelation  function  of 
Hn  n 

the  input  noise 

By  changing  variables  t  =  t  ,  t.-  t  -p  At  ,  Eq.  36  becomes 

k  k  v  k 

E [" C • 3  1  =  (  l/T3)  Z  Z  A(t.  )  A(t.  -p  At  )  p  (p  At  )  At2  = 

L  n  J  ,  k  k  n 

P  k 

=  (  l/Ta  )£  At  pn(p  At)  [i/tE  A(t  )  A(t  -p  At)  At]  (Eq.  37) 
P  n  k 


The  sum  between  the  brackets  in  (Eq.  37)  is  the  time-autocorrelation 

function  R*(pAt)  of  the  sampled  expected  signal.  Although 
s 

depending  on  the  position  of  the  sampling  comb,  its  values  are 

approximately  equal  to  the  corresponding  values  of  the  time- 

autocorrelation  function  R  (p  At)  of  the  non-sampled  signal  (Eq.  25), 

s 
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Thus: 


R^(p  At)  =  0*  [_!-(  Ip  At  |/T)  J  (sin  tt  AF  p  At)/ TT  AF  p  At  (Eq.  38) 

where  |p  At |  <T 

If  the  tapering  factor  l-(|pAt|/T)  is 

E f~C ,a  1  =  (  1/T)E  p  (p  At)  R-  (p  At)  At  =  (  O  2 
L  n  J  n  s  n 

P 

(Eq.  39) 

For  sampling  frequencies  that  are  multiples  of  AF ,  the  sampling 
interval  is 


neglected,  Eq.  37  becomes 


Oga/T)  E(sin  tt  AF  p  At^  AF  p  At)  At 


At  =  l/(m  AF)  (m  integer) 


(Eq.  40) 


Equation  39  then  becomes 

2  mN 

ELCn2  _Kan  as2  T  [sin (  p/m)  tt  /p/m)  TT  J=  (o*  as2  A)/1/™)  L 

p  u  -mN 

^sin  ( p/m)  TT y/  (p/m)  TT  j  ^  Eq  ' 

The  limits  of  p  are  found  from  the  condition  |p  |  At<T-»  |p  |<mT  AF  =mN. 
The  sum  may  be  taken  from  -“>  to  +»,  since  the  tails  contribute 
little  to  its  value  if  N>>1. 
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Since 


L  sin3nx/na  =  £x(tt  -  x)  0  <  x  <  tt  (Ref.  4  )  (Eq.  42) 

n  =  1 


it  follows  that 

CD  /  2  CO 

l/m  £  j_sin  p( Tt/ m)  J  p(Tr/m)  j  =  l/m^  l+2ma/n2Z  '^sinap( n/m)  //p3)] 
P 


=  l/i"^l+|^2ma/TTa]  [in/m(TT-n/m)  ])  =  l/m  [l+(m-l)J  =1  (Eq.  43) 

m  ^  1 


The  power  of  the  noise  at  the  output  of  the 
follows  from  Eq.  41  and  equals 
non-sampling  system  (Eq.  28). 


°s) 


sampling  system  then 
as  in  the  case  of  the 


5 . 3  Signal-to-Noise  Ratio 

Since  neither  the  correlation  peak  nor  the  output  noise  is  affected 
by  the  sampling,  the  output  SNR  of  the  sampled  reference  system  is 
equal  to  the  output  SNR  out  of  the  reference  system  itself.  (See 
also  Appendix  B  of  Ref.  5  ). 


=  N. (SNR) .  (Eq.  44) 

in 


(SNR) 


ref . sampl . 
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NORMALIZED  MODULUS.  NO  PHASE  QUANTIZING 


6  . 


The  first  type  of  quantizing  of  the  received  signal  to  be 
considered  discards  all  amplitude  information  by  normalizing  the 
modulus  of  the  complex  signal,  while  the  phase  information  is 
retained  completely. 


Technically  this  type  of  quantizing  occurs  when  a  high-frequency, 
narrow-band  signal  is  clipped.  After  the  clipping  it  is  possible  — 
by  filtering  and  analogue  quadrature  demodulation  —  to  extract 
the  phase  of  the  modulating  low-frequency  function  with  complete  loss 
of  amplitude  information.  The  accuracy  of  the  phase  approximation 
is  a  function  of  the  relative  bandwidth  Q  ~  F/AF  of  the  high- 
frequency  signal.  The  results  of  this  chapter  are  valid  for  a 
technical  system  in  which  F/AF>>1. 


If  the  modulus  of  the  received  complex  signal  is  normalized  to 

S  =  c  J2 ,  we  may  write  (Fig.  4) 
s 


£  (t) 

norm 

with 


jarc  £(t) 

b  e 

£(t)  =  aS  e^s^^  +v(t) 
v(t)  =  N(t)  e^n^)  =  input 


noise 


(Eq.  45) 
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6 .  1  Signal  Correlation  Peak 


The  mean  value  of  the  correlation  peak  is 


e[y(o) ]  =  e[i/t 


&  (t) 

norm 


=  S3/T  jB[eJarcUt)]  e-^(t)  dt 


(Eq.  46) 


Since  the  noise  vector  v(t)  has  a  uniform  phase  distribution, 

the  phase  of  E  E,  (t)  will  be  equal  to  the  phase  cp  (t)  of  the 

L  norm  J  s 

expected  signal.  By  pairing  positive  and  negative  values  of  the 

phase  cp  .  around  cp  (t),  it  follows  for  a  given  value  of  N(t) 
noise  s 

(see  Fig.  5)  that 


E|-ejarc5(t)|N(t)-|  .  edcps(t)  E[cos(arc,r(t))  |M(t)] 
,  .  __  aS  +  N  cos  cp 


i  r 


Cp  =0  J 
n 


.  J op„ ( t )  /aS  \ 

dcpn=e  s  1  irr) 


a3S2  +  2aSNcos  cp  +N3 
n 


(Eq.  47) 


where 


TT 


T(~—J 

\  N  7 


JL 

TT 


aS/  +  cos  cp 
N  ^n 


dcp 


(Eq.  48) 


0  y(aS/N)5+  1  +  2^aS//N)cos  cp 
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Equation  48  can  be  solved  using  complete  elliptic  integrals,  but 

aS  aS 

the  limiting  cases  of  ~  «1  and  ~  >>1  are  readily  found 
1  N  N 

directly: 


.  ^  aS 

lf  TT 


»1 


(Eq.  49) 


To  find  the  average  E^e^arC  J,  the  integral  in  Eq.  48  has  to  be 

averaged  over  N(t) ,  ranging  from  0  to  »  with  a  Rayleigh 
distribution.  Thus 


E[eJarc  5<t)j  _ 


(Eq.  50) 


Substitution  in  Eq.  46  yields  the  mean  value  of  the  complex 
correlation  peak  of  the  modulus-normalizing  system.  For  the  limiting 
cases : 


a  a 


n 


aa 

if  — -  «1 

a 

n 


aa 

if  - -  »1 

a 

n 


(Eq.  51) 


If  the  system  is  sampled,  the  same  results  hold. 
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6 . 2  Noise 


In  the  absence  of  signal,  the  correlator  output  at  a  given  instant  t, 
equals 


/.  n  _  1  r  o  “JCP  (fc)  j.  S  P  jcp  (t)  -jcp  (t) 

v  (t,  =  “  £,  (.  t )  S  e  s  dt=  ~  I  e  n  e  s  dt 

Tn  '  T  J  norm  T  J 

(Eq.  52) 

This  integral  is  a  complex  random  variable  with  zero  mean,  since 
the  mean  of  the  normalized  noise  vector  e^n^^  equals  zero. 

The  output  power  is  equal  to  the  variance  of  the  complex  variable 
(Eq.  52) 


r  j  *  2  ~l 

r  *i  s4  f 

T  r  j[<Pn(t)-*n(9)]- 

-j[cps(t)-tps(  e)  ] 

Lki2} 

LYn  Yn  ]“  a  Ji 

J  EL* 

]  * 

dt  d0 


s_ 

T 


dT  Ej-eJtCP«(t)'0Pr1(t-T)  j'  r,  je-j[vs(t)-Vs(t-T)] 


dtj  (Eq.  53) 


with  T  =  t  -  0 


The  integral  between  brackets  in  the  last  expression  of  Eq.  53  is  the 
conjugate  complex  time-autocorrelation  function  of  the  expected 
complex  signal  This  is  a  sin  tt  AF  t/tt  AF  t  function 

with  zero  phase,  since  the  energy  spectrum  of  the  expected  signal  is 

\ 

a  realj f lat  band  of  ±^AF  around  zero  frequency.  See  also  Ref.  5 
Appendix  A.  The  tapering  factor  l-( |t|/T)  due  to  the  limited 
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time-duration  of  the  signal  will  be  neglected,  as  in  the  foregoing 
chapters.  Also,  the  integration  limits  -T  AF  to  +T  AF  will 
be  extended  to  -«  to  +»,  since  the  tails  of  the  integration 
contribute  little  to  the  solution  of  Eq.  53  if  T  AF  =  N>>1.  Thus 


4  n 


E[>nl  ]  H  J 


d(  AF  t) 


sin  n  AF  T 

n  AF  t 


P  (t) 

rn  norm 


(Eq.  54) 


|(  T )  =  E[( 


3'[cpn(t)  -  cp„(t  -  T) 


is  the  auto-correlation 


pn  norm 

function  of  complex  Gaussian  noise  normalized  to  unit  modulus*  The 
correlation  function  may  be  written  as 


n  norm 


r  j(ept-epa) 

(t)  =  B[o  ]  - 


2tt  2tt  . , 

p  J  ( cp ,  -cpa  ) 

e 


Cp ,  =0  cpa=0 


2tt 

r> 


2TT-CPa 


p( Cp ,  ,Cpa)  dcp,  dcpa 


(Eq.  55) 


dcpa 


eJ  ^  p[f  +  Cpa,  Cpa  ]  df  , 


cps=0  t=Cpa 


where  t  =  cp,  -  cp2 
and  (Ref.  3  p.  164): 


>LCp,  ,CPa]: 


°n-p„a'T> 


4rra  a  4 

n 


t 


(l-P3)2  +  p(rt-arc  cos  p) 


(1-03) 


ax3/2 


]  =  PU)  (Eq.  56) 
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in  which 


Pn(T) 


cos  (Qp,-qp2) 


sin  tt  AF  t 

TT  AF  t 


cos  l|f 


p  (  T)  _  q  a - —  =  autocorrelation  function  of  noise 

n  n  TT  AF  T 

components 

The  function  p j^cp ,  ,  cpa  j  is  the  joint  probability  density  function 
of  qp(  ,  qpa  and  is  an  even  function  p(\(r)  of  cp ,  —  Qpa  ~  ♦  >  as 
follows  from  (Eq.  56).  If  t  =  0,  Eq.  56  should  reduce  to 


PoU) 


(Eq.  57) 


since  the  phase  difference  between  the  noise  vectors  is  then  zero 
with  probability  one. 


The  last  integrand  in  Eq.  55  is  periodic  in  ^ ,  so  that  the 
integration  may  be  taken  over  the  interval  ( — tt,  +tt)  : 

TT 

p  (t)  =  2tt  f  (cos  i/  +  j  sin  if)  p (if)  df  (Eq.  58) 

rn  norm  J 

-TT 


The  term  with  sin  ij)  will  give  zero,  since  p(^)  is  an  even 
function  of  i|r .  Thus: 


P 


n  norm 


(  t) 


TT 

=  2tt  cos  i|i  p(i|t)  di|r 

C 

-TT 


(Eq.  59) 
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It  follows  from  Eq.  57  that  for  t  =  0 


P  (0) 

n  norm 


=  J  6(f)  d*  -  1  , 

-tt 


(Eq.  60) 


as  was  to  be  expected. 


The  noise-output  power  becomes,  from  Eqs.  54,  56,  &  59: 


r  3 1  s4  r 

L  I  Yfi  I  J  =  N~  dx  sinC  x  (1  -  sine2  x) 


1 

2n  J 

-TT 


M  r~ 

(l-sinc2  xco^  f)2+sinc  xcos  ^  TT-arc  cos(sinc  x  cos  i|r ) 
di|f  cos  iji  — - - - - 


(l  -  sine8  x  cos2  i|i ) 


3/2 


with  (  sine  x  = 


(Eq.  61) 


sin  ttx  ,  _  _ 

- 1 -  (Ref.  2  ) 

TTX 


x  =  AF  t 

The  solution  of  Eq.  6l  has  not  been  found,  so  that  the  SNR  at  the 
output  of  the  modulus-normalizing  system  cannot  be  given. 

If  the  system  is  sampled  with  a  sampling  interval 


At 


m  AF 


(m  integer) 
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it  can  be  shown,  by  following  the  same  reasoning  as  in  Chapter  5> 
that  Eq.  54  takes  the  form 


n  i2i 

.§1 

mN 

P 

sme  - 

P  (  P  ) 

UYP  j 

N 

m 

L 

-mN 

m 

n  norm  m  AF  J 

(Eq.  62) 


in  which 


if  p  =  0  p  (0)  =  1 

rn  norm 


(Eq.  63) 


if  p  7^  0  p 


n  norm  \m 


V — )  =  ( 1-  sine2  £  )  • 
'm  AF  /  \  m  J 


TT 


(*- 


2d  2  _l e-i —  £  cos  f  JjT  -arc  cos  ^sine  ^  cos 


sine  2£  cos2  to )  2  +sinc  i:  co 
m  m 


£ 

m 


*  2tt  •  w 

dijr  cos  i|r  " 

TT 

(l-  sine2^  cos2  *  )3/2 

From 

m  =  1 

( sampling 

with  frequency  AF),  Eq * 

62  becomes 

E[l 

-< 

II 

N 

^  sin  p  tt 

p=-N  p  TT 

P  i-2-) 

n  norm  \^p  / 

(Eq.  64) 

From  Eq.  63  it  follows  that 


f— 

^n  normVAF 


so  that 


1  if  p  =  0 
0  if  p  ^  0 


(Eq.  65) 


r  3  -1  S* 

E[lvJ  ]-  V- 


(sampling  freq.  =  AF) 


(Eq.  66) 
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6 . 3  Signal-to-Noise  Ratio 


The  output  SNR  of  the  modulus-normalizing  system,  sampled  with 
frequency  AF,  becomes  (Eqs.  51  &  66) 


N *  j  (-^“)  =  0.785  •  N  ‘  (SNR) 
n 


a  a 


if  - 

in  a 


<<  1 


n 


(SNR) 


norm 


(Eq.  67) 


ao 


>>1 


n 


It  is  seen  that  such  a  system  is  about  1  dB  lower  than  the  reference 
system  for  low-input  SNR. 

For  high  input  SNR,  the  output  SNR  of  the  modulus-normalizing 
system  saturates  at  a  value  N. 
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7. 


90° -PHASE  QUANTIZING 


If  four  standard  vectors  making  angles  of  90°  are  used  to  quantize 
the  complex  input  signal,  the  output  properties  will  be  independent 
of  the  orientation  of  the  standard  vector  star.  This  is  a 
consequence  of  the  assumed  uniform  phase  distribution  of  the  signals 
and  the  complete  phase  description  of  the  expected  signal. 

We  will  determine  the  output  SNR  of  a  system  with  the  vector  star 
orientated  at  45°  with  respect  to  the  axes  of  the  complex  plane. 

(Fig.  6). 


Fig.  6 

The  real  and  imaginary  parts  of  the  complex  input  signal  can  then 

only  take  the  values  +ff  and  -a  i.e.  they  are  infinitely  clipped. 

s  s 
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The  following  notation  will  be  used  for  the  clipped  components: 


+c  if  a(t)  >  o 

s 

5t(t)  =  (Eq.  68) 

-c  if  a(t)  <  0 

s 

7 . 1  Signal  Correlation  Peak 

If  cp„  =  0 ,  the  cosine  component  of  the  input  is 

a(t)  =  a  A(t)  +  n(t) 

The  mean  value  of  the  correlation  peak  is 

E[CM0)]  =  “  f  E^a(t)J  A(t)  dt  (Eq.  69) 

To  determine  E^a(t)j,  we  consider  two  cases,  A(t)>0  and 
A(t)<0,  it  follows  then  that 

A(t)>0  -»  Ej^a(t)  1  =  (  +  as)Probj^n(t)  >-a  A(t)  J  +  (-ag)Probj^n(t)<-a  A(t)j 

A(t)<0  -»  Ej^a  ( t )  J  =  (+0  )  Probj^n(  t )  >|a  A(  t )  |  J  +(  )  Prob^n(  t  )<  |a  A(t)|j 

(Eq.  70) 
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Fig. 7 


It  follows  from  Eq.  70  and  the  Gaussian  amplitude  distribution 
of  the  noise  (Fig.  7)  that  the  probability  that  the  mean  value  of 
the  clipped  signal  at  a  given  instant  has  the  same  sign  as  that  of 

i 

the  expected  signal  A(t)  is 


lj>  lA(t) 
L  a 

n 


a|  A(t)|  -2q  2 

n 


a  lA(t)  1 


on  r*r  j 


e  dn  =  ,/2/tt 


-a|  A(t)[ 


n 


-K; 


d^ 


(Eq.  71) 


The  mean  value  of  the  clipped  signal  thus  becomes: 


(Eq.  72) 


Two  limiting  cases  will  be  considered:  low  and  high  input  SNR. 
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that 


tp  a  lA(t)  1 

a*  ir  «  1 ,  it  follows  from  a  power  series  expansion 

n 


pjALLU  J  =  j  2/n 

n 


lA(t) 


(Eq.  73) 


n 


Thus,  from  Eq.  72, 


{«<*)]  -  °s  v  2/tt  a|^(t)l  -  /27n  Aft)  (Eq.  74) 

J  | Aft ) |  s 


a  a 

n  n 


So  that  the  correlation  peak  (Eq.  69)  becomes 


Lca(0)]  =  i  |As(t)dt  = 


(Eq.  75) 


b.  If 


a|A(t)  | 


>>1,  then  I 


aJA(.t)  ll  1. 


n 


u  a 


n 


Consequently  E^a(t)j 


=  A(t) 
I A  ( t )  I 


a  ,  i.e.  the  signal-plus-noise 


has,  with  probability  one,  the  same  sign  as  the  signal  only,  as  is 
obvious.  The  mean  value  of  the  correlation  peak  becomes  in  this 
case,  from  Eq.  69: 


E[CS(0)]  -  T 


’a  TmoI  A(t)  dt  =  ~r 


|  A ( t )  I  dt  (Eq.  76) 
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The  last  expression  involves  calculation  of  the  average  absolute 
value  of  the  cosine  component  A(t)  of  the  expected  signal.  Because 
of  the  uniform  phase  distribution  and  the  constant  amplitude, 
the  average  absolute  value  is  equal  to  that  of  a  sinusoid  with  the 
same  amplitude.  Thus 


1 

T 


"*  2 

|A(t)ldt=~a  a/2 
1  1  tt  s 


(Eq.  77) 


where  a  *J2  =  amplitude  of  the  signal  component. 

S 

So  that 


nT1 

> 

o 

ll 

0, 

L  S  J 

TT  S 

(Eq.  78) 


If  the  input  signal  is  not  only  quantized,  but  also  sampled,  the 
same  results  are  valid.  See  also  Chapter  5» 


7.2  Noise 


In  the  absence  of  signal,  the  correlator  output  at  a  given  instant 
is 


n(t)  A(t)  dt 


(Eq.  79) 


This  equation  has  the  same  form  as  Eq.  21.  The  variance  of  the 
random  variable  (Eq.  79)  can  be  determined  by  following  the  same 
reasoning  as  in  Chapter  4.2. 


44 


Taking  into  account  that  the  autocorrelation  function  of  the 
clipped  Gaussian  noise  is 


Pfi<T) 


2 

TT 


■  P  (  T) 


~  a  3  arc  sin  [~— - "| 

8  L  o=  J 

n 


2  3 

~  a  arc  sin 
tt  s 


sin  tt  AF  T  1 
L  TT  AF  T  J  * 


(Eq.  80) 


the  expectation  of  the  squared  correlator  output  for  noise  only 
becomes  (Eqs.  22  to  27): 


T  AF 


|~r  al  1  _  4 

sin  TT  AF  t 

2  •  r 

arc  sin 

TT  L 

sin  TT  AF 

L  n  j  T  AF  s  v 

TT  AF  T 

3 

1 _ 

-T  AF 


(Eq.  81) 


The  integral  in  the  expression  above  has  been  calculated  to  be 
0.7898  for  the  limits  -»  to  +».  The  noise  power  at  the 
correlator  output  becomes 


cT1 

> 

to 

I! 

0.7898  ^  o4 

L  n  J 

N  s 

(Eq.  82) 


If  the  system  is  sampled  with  a  sampling  interval  of 


At  = 

X 

m  AF 

(m 

integer) , 

Eq.  81 

takes 

the 

form 

r?  1  fy  •  3  I 

-  1  n  4 

1 

mN 

z 

p=-mN 

P 

sin  ^  17 
m 

2  •  1 

•  “  arc  sin 
TT  1 

sin 

r 

E  17 
m 

i-k 

ELcs  j 

N  °s 

m 

^  17 
m 

L 

E 

17 

J  k  N 

m 


with  N  =  T . AF 


(Eq.  83) 
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The  results  in  Table  1  were  found  for  the  factor  k  in  Eq.  83 
when  N>>1: 


m 

m 

1  S  =k 

m 

1  2 

k  TT 

1  H 
k  4 

1  8 
k  tt3 

2 

k 

1 

1.000 

0.6366 

0.7855 

0.8105 

1.0000 

2 

0.8402 

0.7577 

0.9348 

O.9648 

1. 1902 

3 

0.8117 

0.7844 

0.9676 

0.9986 

1.2320 

4 

0.8021 

0.7937 

0.9794 

1.0106 

1.2467 

8 

0.7929 

O.8O3O 

0.9906 

1.0223 

1.2612 

16 

0.7906 

0.8052 

0.9935 

1.0252 

1.2648 

32 

0.7901 

0.8058 

0.9940 

1.0260 

1.2657 

64 

0.7899 

0.8060 

0.9942 

1.0262 

1.2660 

128 

0.7898 

0.8062 

0.9945 

1.0264 

1.2663 

Noise  figure 

TABLE  1 

k  and  output  SNR  for  low  and  high  input 

SNR  as  a 

function  of  the  sampling  rate 

It  is  seen  from  Table  1  that  the  noise  figure 


sin  ^  tt  sin  ^  tt 

m  2  m 

-  .  —  arc  sin  f—  ■  -  —  j  (Eq.  84) 

TT  L  p  J 

p  “  rr 

—  tt  m 

m 

converges  rapidly  to  an  asymptotic  value  of  about  0.7898,  as  given 
in  Eq,  82  for  the  solution  of  the  integral  in  Eq.  81. 
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7 . 3  Slgnal-to-Noise  Ratio 


The  output  SNR  of  the  90°-phase  quantizing  correlator  is  found 
from  Eqs.  75.  78,  &  83  and  Table  1. 


For  low  and  high  input  SNR: 


(SNR) 


out 


f 


if  (SNR) .  « 1 


(Eq.  85) 


if  (SNR).  »1 
m 


It  follows  from  Table  1  that,  for  normal  sampling  (AF)  resp, 
infinite  sampling  the  losses  as  compared  to  the  reference  system  vary 
from  2  dB  to  1  dB  for  low  input  SNR.  For  high  input  SNR  the 
performance  under  normal  sampling  is  1  dB  less  than  that  of  a 
correlator  with  no  phase  quantizing  (Chapter  6.3). 

The  behaviour  of  the  non-sampling  system  cannot  be  compared  to 
the  non-sampling,  non-phase  quantizing  system  of  Chapter  6,  since 
the  solution  of  Eq ,  6l  has  not  been  found. 
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8. 


-  -PHASE  QUANTIZING 
m 


If  the  complex  input  signal  is  regularly  phase-quantized  to  angles 
smaller  than  90° ,  it  is  for  some  cases  possible  to  determine  the 
effects  of  this  finer  quantizing  on  the  signal  correlation  peak 
and  on  the  noise , 

8 . 1  Signal  Correlation  Peak 

Independent  of  the  sampling  frequency,  the  mean  signal  correlation 
peak  can  be  determined  directly  if  the  input  SNR  is  high.  In  that 
case  the  received  complex  signal  will  fall,  with  probability  one, 
in  the  sector  belonging  to  the  expected  complex  signal  (Fig.  8). 


\ 


Real 


Fig.  8 


Fig.  9 
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The  phase  difference  between  the  quantized  received  vector  and  the 

non-quantized  expected  signal  vector  ranges  therefore  from  — 

TT  m 

to  +£  “  .  (See  Fig.  9).  On  the  average,  the  length  in  the 
direction  of  the  correlation  peak  (phase  angle  zero,  since  it  is 
assumed  that  cp„  =  0)  becomes  smaller  with  a  factor; 


TT 

m 


JL 

2m 

cos  cp  dcp 

TT 

2m 


sin 


rr 

2m 


JL 

2m 


(Eq.  86) 


The  modulus-normalizing  system  of  Chapter  6  should  be  taken  as 
a  reference.  It  was  shown  there  (Eq.  51)  that  for  high  input  SNR 
the  signal  correlation  peak  saturates  at  a  value  2o  2  . 


Thus  the  mean  value  of  the  correlation  peak  is 


sin  77 

2m 

2a  2 

s 

TT 

2m 

(Eq.  87) 


As  an  example  we  consider  the  case  m  =  2,  i.e.  90°-phase  quantizing. 


Then 


2o 2  -  2^2 

s  TT 


(Eq.  88) 
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Compare  this  result  with  Eq.  78,  which  is  half  its  value  because 
it  gives  the  contribution  of  only  one  component  to  the  correlation 
peak,  (a  A) . 


8 . 2  Noise 

If  the  sampling  frequency  is  AF,  the  complex  noise  samples 

are  uncorrelated.  Independent  of  the  phase-quantizing,  the  output 

of  the  correlator  will  be  the  sum  of  N  vectors  with  equal  lengths 

and  random  orientation.  For  the  limiting  case  of  infinitely 

fine  phase  quantizing,  the  result  has  already  been  found  in  Eq.  66.: 


4 


s 


(Eq.  89) 


N 


8 . 3  Signal-to-Noise  Ratio 

For  a  phase-quantizing  correlator  sampling  with  sampling  frequency  AF, 
the  output  signal  SNR  for  high  input  SNR  is 


sin  —II—  2 


(SNR) 


out 


(Eq.  90) 


2m 
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9 


COMBINATION  OF  RESULTS 


The  results  found  in  the  foregoing  chapters  will  now  be  combined 
and  shown  graphically  as  functions  of  the  degree  of  phase-quantizing, 
the  sampling  interval,  and  the  input  signal-to-noise  ratio. 

9 . 1  Signal  Correlation  Peak 

The  power  of  the  signal-output  peak  is  represented  in  two 
three-dimensional  drawings:  Fig.  10  for  low  input  SNR,  Fig.  11  for 

4 

high  input  SNR.  A  scaling  factor  N/4d  is  used  to  normalize 

s 

the  results  of  Chapters  6,  7,  and  8. 
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For  low  input  SNR,  the  signal  peak  power  is  independent  of  the 
sampling  interval,  varying  from  rr/4  N.(SNR)in  for  no-phase 
quantizing,  to  2/tt  N.(SNR)^n  for  90°-phase  quantizing.  The 
intermediate  values  probably  follow  the  function 

rr/4  [sin(  tt/ 2m)  /  (lT/2m)]2  N.(SNR)±n 


Fig.  11 
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For  high  input  SNR,  the  signal  peak  power  saturates  to  a  value 
N  in  the  case  of  no-phase  quantizing  (only  modulus  normalizing). 
With  increasing  phase  quantizing,  the  signal  level  goes  down 
according  to  a  ^sin(ir/2m)  j  (Tr/2m)J  -  law,  where  ir/m  is  the 
step  of  the  phase  quantizing.  The  output  peak  power  does  not 

4 

depend  on  the  sampling  frequency.  The  scaling  factor  is  40s  • 

9 . 2  Noise 

The  output  noise  power  for  a  modulus  normalizing,  phase  quantizing, 
and  sampling  complex  correlator  is  given  in  Fig.  12. 

The  values  found  in  Chapters  6,  7,  and  8  are  multiplied  by  the 

/  4 

scaling  factor  N/4CT  to  make  the  results  comparable  to  the 

s 


signal  values  of  function 


Outpub  noise  level 


0 


±_ 

AF 


Samplt  nty  i 


interval  af 


Phase  -ft 

quantizing  — 


Fig.  12 
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In  Chapter  8  it  has  been  shown  that  when  sampling  with 
frequency  AF ,  the  noise  output  power  is  independent  of  the 
phase  quantizing.  Although  not  proved,  it  is  assumed  that  for 
other  sampling  frequencies  the  output-noise  power  is  also 
independent  of  the  phase  quantizing.  This  means  that  for  no-phase 
quantizing  (only  modulus  normalizing)  the  noise-output  power  as 
a  function  of  the  sampling  rate  is  that  given  in  Table  1  (dotted 
line  in  Fig.  12) . 

9 . 3  Signal-to-Noise  Ratio 

The  output  signal-to-noise  ratios  for  low  and  high  input  SNR 
(Figs.  13  and  14)  follow  directly  from  the  above  paragraphs. 


Fig.  13 
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If  we  accept  the  assumption  that  the  output-noise  power  is 
independent  of  the  phase  quantizing,  then  the  output  SNR  for 
modulus  normalizing  only  is  about  equal  to  the  output  SNR  of  the 
reference  system  for  low  input  SNR.  Both  quantizing  and  sampling 
of  the  complex  input  signal  reduce  the  output  SNR  by  about  1  dB. 

If  the  input  signal  is  not  only  quantized  by  90° ,  but  also  sampled 
with  sampling  frequency  AF,  the  total  loss  is  about  2  dB. 


For  high  input  SNR,  the  output  SNR  saturates  at  a  maximum  value 
of  1.266  N  (Table  l) . 
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In  Fig.  15  the  combined  influences  of  modulus  normalizing,  phase 
quantizing,  and  sampling  are  given  as  functions  of  the  input  SNR. 
The  curves  are  interpolated  between  low  and  high  input  SNR. 
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It  follows  from  Fig.  15  that  normalizing  only  the  modulus 
causes  saturation  of  the  input  SNR  for  high  input  SNR. 

Sampling  and  no-phase  quantizing,  or  90°-phase  quantizing  and  no 
sampling,  makes  the  system  about  1  dB  less  than  the  modulus 
normalizing  system. 

Both  90°-phase  quantizing  and  sampling  cause  a  total  loss  of 
2  dB. 

It  has  been  shown  in  Ref.  5  that  90°-phase  quantizing  of  the 
expected  signal  causes  an  additional  loss  of  1  dB.  This  would 
raise  the  losses  to  a  total  of  3  dB  as  compared  with  those  of 
the  modulus  normalizing  system  without  phase-quantizing  of  the 
expected  signal,  and  without  sampling. 
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SCNR  Turkey  1 
SCNR  U.K.  1 
SCNR  U.S.  1 

NLO  to  SACLANTCEN 

NLO  France  1 
NLO  Italy  1 
NLO  Portugal  1 
NLO  U.S.  1 
NLO  U.K.  1 

NLR  to  SACLANT 

NLR  Belgium  1 
NLR  Canada  1 
NLR  Denmark  1 
MM  France  1 
NLR  Germany  1 
NLR  Greece  1 
NLR  Italy  1 
NLR  Norway  1 
NLR  Portugal  1 
NLR  Turkey  1 
NLR  U.K.  1 


